The Journal of Physical Chemistry Letters

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump probe experiments for investigating electronic and nuclear changes during light-induced reactions. On time scales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer1. However, all ultra-fast TR-SFX studies to date have employed such high pump laser energies that several photons were nominally absorbed per chromophore2-14. As multiphoton absorption may force the protein response into nonphysiological pathways, it is of great concern15 whether this experimental approach16 allows valid inferences to be drawn vis-a-vis biologically relevant single-photon-induced reactions17. Here we describe ultrafast pump-probe SFX experiments on photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe-CO bond distance (predicted by recent quantum wavepacket dynamics15) are seen to depend strongly on pump laser energy. Our results confirm both the feasibility and necessity of performing TR-SFX pump probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing design and interpretation of ultrafast TR-SFX pump probe experiments16 such that biologically relevant insight emerges.

The DNA sequences available in the prebiotic era were the genomic building blocks of the first life forms on Earth and have therefore been a matter of intense debate.1,2 On the surface of the Early Earth, ultraviolet (UV) light is a key energy source3, which is known to damage nucleic acids4. However, a systematic study of the sequence selectivity upon UV exposure under Early Earth conditions is still missing. In this work, we quantify the UV stability of all possible canonical DNA sequences and derive information on codon appearance under UV irradiation as selection pressure. We irradiate a model system of random 8mers at 266 nm and determine its UV stability via next-generation sequencing. As a result, we obtain the formation rates of the dominant dimer lesions as a function of their neighboring sequences and find a strong sequence selectivity. On the basis of our experimental results, we simulate the photodamage of short proto-genomes of 150 bases length by a Monte Carlo approach. Our results strongly argue for UV compatibility of early life and allow the ranking of codon evolutionary models with respect to their UV resistance.

Chlorophylls (Chls) are known for fast, sub-picosecond internal conversion (IC) from ultraviolet/blue absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates, such as photosynthetic pigment-protein complexes. We show that according to a simple Forster resonance energy transport (FRET) scheme, unmodulated B-B EET likely poses an existential threat, in particular the photochemical reaction centers (RCs). This insight leads to so-far undescribed roles for carotenoids (Crts, cf. previous article in this series) and Chl b (this article) of possibly primary importance. It is demonstrated how pigments in a photosynthetic antenna pigment-protein complex (CP29) undergo FRET. Here, the focus is on the role of Chl b for EET in the Q and B bands. Further, the initial excited pigment distribution in the B band is computed for relevant solar irradiation and wavelength-centered laser pulses. It is found that both accessory pigment classes compete efficiently with Chl a absorption in the B band, leaving only 40% of B band excitations for Chl a. B state population is preferentially relocated to Chl b after excitation of any Chls, due to a near-perfect match of Chl b B band absorption with Chl a B state emission spectra. This results in an efficient depletion of the Chl a population (0.66 per IC/EET step, as compared to 0.21 in a Chl a-only system). Since Chl b only occurs in the peripheral antenna complexes, and RCs contain only Chl a, this would automatically trap potentially dangerous B state population distantly from the RCs.

Chlorophylls (Chls) are known for fast, sub-picosecond internal conversion (IC) from ultraviolet/blue absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates, such as photosynthetic pigment-protein complexes. We show that according to a simple Forster resonance energy transport (FRET) scheme, unmodulated B-B EET likely poses an existential threat, in particular the photochemical reaction centers (RCs). This insight leads to so-far undescribed roles for carotenoids (Crts, this article) and Chl b (next article in this series) of possibly primary importance. Here we show that B [->] Q IC is assisted by the symmetry-allowed Crt state (S2) by using the plant antenna complex CP29 as a model: The sequence is B [->] S2 (Crt, unrelaxed) [->]S2 (Crt, relaxed) [->] Q. This sequence has the advantage of preventing ~ 39% of Chl-Chl B-B EET, since the Crt S2 state is a highly efficient FRET acceptor. The likelihood of CP29 to forward potentially harmful B excitations towards the photosynthetic reaction center (RC) is thus reduced. In contrast to the B band of Chls, most Crt energy donation is energetically located near the Q band, which allows for 74/80% backdonation (from lutein/violaxanthin) to Chls. Neoxanthin, on the other hand, likely donates in the B band region of Chl b, with 76% efficiency. The latter is discussed in more detail in the next article in this series. Crts thus do not only act in their currently proposed photoprotective roles, but also as a crucial building block for any system that could otherwise deliver harmful "blue" excitations to the RCs.

Retinal is a chromophore covalently bound to various photoreceptors. Its photo-induced isomerization triggers a series of structural changes named photocycle, leading to diverse biological functions. Despite tremendous advances in structural biology and artificial intelligence-driven structure prediction, it remains challenging to analyze all photocyclic intermediates. Here, we present an optimized computational approach to calculate RSBH+ isomerization and its induced structural changes based on a classical molecular mechanics approach using quantum mechanically improved retinal force field parameters. Isomerization is induced by an excited state restraint which is subsequently relaxed to allow the return to the electronic ground state. We applied this approach to the key protein of optogenetics, Channelrhodopsin-2 from Chlamydomonas reinhardtii (CrChR2). Besides the reformation of the alltrans/CN-anti ground state, we observed the production of a mixture of two isomeric states 13-cis/CN- anti and 13-cis/CN-syn. These findings agree with the previously found branched photocycle model based on experimental data. Our calculations show an asymmetric potential energy landscape of the excited state leading to a corresponding isomerization state distribution. Unlike earlier publications, our procedure describes the retinal photoisomerization on the natural timescale of 500 fs. As our newly derived retinal force field parameter set precisely relies on quantum biological knowledge, it assists to improve the refinement of experimental structure biological data. Our readily customizable strategy provides mechanistic insights at high spatio-temporal resolution, which permits accurate structural predictions of early photocycle intermediates. These insights will stimulate the rational design of optogenetic tools thus providing improved diagnostic and therapeutic approaches for neuronal and other diseases. HighlightsO_LIuniversal method to study molecular mechanism of optogenetic tools C_LIO_LIretinal photo-isomerization calculation in real time C_LIO_LIprediction of branched photo cycle agrees with experimental IR spectroscopic results C_LIO_LIdetected asymmetric excited state potential energy landscape C_LIO_LIassists to improve structural model refinement of retinal proteins C_LI Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=116 HEIGHT=200 SRC="FIGDIR/small/707937v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@71f7fdorg.highwire.dtl.DTLVardef@503482org.highwire.dtl.DTLVardef@1a77120org.highwire.dtl.DTLVardef@1f410a0_HPS_FORMAT_FIGEXP M_FIG C_FIG

The accumulation of insoluble protein aggregates containing amyloid fibrils has been observed in many different human protein misfolding diseases1,2, and their pathological features have been recapitulated in diverse model systems3. In vitro kinetic studies have provided a quantitative understanding of how the fundamental molecular level processes of nucleation and growth lead to amyloid formation4. However, it is not yet clear to what extent these basic biophysical processes translate to amyloid formation in vivo, given the complexity of the cellular and organismal environment. Here we show that the aggregation of a fluorescently tagged polyglutamine (polyQ) protein into {micro}m-sized inclusions in the muscle tissue of living C. elegans can be quantitatively described by a molecular model where stochastic nucleation occurs independently in each cell, followed by rapid aggregate growth. Global fitting of the image-based aggregation kinetics reveals a nucleation rate corresponding to 0.01 h-1 per cell at 1 mM intracellular protein concentration, and shows that the intrinsic stochasticity of nucleation accounts for a significant fraction of the observed animal-to-animal variation. Our results are consistent with observations for the aggregation of polyQ proteins in vitro5 and in cell culture6, and highlight how nucleation events control the overall progression of aggregation in the organism through the spatial confinement into individual cells. The key finding that the biophysical principles associated with protein aggregation in small volumes remain the governing factors, even in the complex environment of a living organism, will be critical for the interpretation of in vivo data from a wide range of protein aggregation diseases.

Photosystem I (PSI) contributes to light-conversion reactions; however, its oligomerization state is variable among photosynthetic organisms. Herein we present a 3.8-[A] resolution cryo-electron microscopic structure of tetrameric PSI isolated from a glaucophyte alga Cyanophora paradoxa. The PSI tetramer is organized in a dimer of dimers form with a C2 symmetry. Different from cyanobacterial PSI tetramer, two of the four monomers are rotated around 90{degrees}, resulting in a totally different pattern of monomer-monomer interactions. Excitation-energy transfer among chlorophylls differs significantly between Cyanophora and cyanobacterial PSI tetramers. These structural and spectroscopic features reveal characteristic interactions and energy transfer in the Cyanophora PSI tetramer, thus offering an attractive idea for the changes of PSI from prokaryotes to eukaryotes.

Quinone-based electron bifurcation (EB) catalyzed by cytochrome bc1 (cytbc1) plays a critical role in maximizing efficiency of biological energy conversion. The canonical EB model (CEB), grounded in equilibrium redox potentials, dictates the order of EB steps with initial endergonic reduction of high-potential iron-sulfur cluster (2Fe2S) by quinol followed by exergonic reduction of low-potential heme bL (bL) by semiquinone (SQ). However, this concept falls short in explaining several experimental observations, including intermediate semiquinone spin-coupled to 2Fe2S (SQ-2Fe2Sred) and the absence of short-circuiting. Presented here DFT calculations on large cluster models of cytbc1, encompassing both 2Fe2S and bL, identified location of donor (HOMO) and acceptor (LUMO) orbitals along with the previously not considered microstates to reveal that EB is an emergent property of an integrated system of redox cofactors where transient charge separations dynamically modulate electron affinities. In this system, electron transfer initiates preferentially toward bL, indicating a departure from the conventional sequence proposed by CEB. Based on this finding, we introduce an EMBER (EMergent BL-first Electron Routing) model of EB and demonstrate that its assumptions are supported by electron paramagnetic resonance spectroscopy data. Unlike CEB, EMBER proposes a relatively flat energy profile for EB that accommodates stable SQ-2Fe2Sred and explains suppression of short-circuits without additional assumptions. It highlights the importance of state-dependent electrostatic interactions in shaping electron transfer pathways. In general, the concept of emergence inherent to EMBER offers a mechanistic framework applicable to a broad range of multi-cofactor redox enzymes beyond cytbc1.

The structural dynamics of a molecule are determined by the underlying potential energy landscape. Conical intersections are funnels connecting otherwise separate energy surfaces. Posited almost a century ago 1, conical intersections remain the subject of intense scientific investigation 2-4. In biology, they play a pivotal role in vision, photosynthesis, and DNA stability 5,6. In ultrafast radiationless de-excitation 1,7, they are vital to ameliorating photon-induced damage. In chemistry, they tightly couple the normally separable nuclear and electronic degrees of freedom, precluding the Born-Oppenheimer approximation 8. In physics, they manifest a Berry phase, giving rise to destructive interference between clockwise and anti-clockwise trajectories around the conical intersection 9. Accurate theoretical methods for examining conical intersections are at present limited to small molecules. Experimental investigations are challenged by the required time resolution and sensitivity. Current structure-dynamical understanding of conical intersections is thus limited to simple molecules with around 10 atoms, on timescales of about 100 fs or longer 10. Spectroscopy can achieve better time resolution, but provides only indirect structural information. Here, we present single-femtosecond, atomic-resolution movies of a 2,000-atom protein passing through a conical intersection. These movies, extracted from experimental data by geometric machine learning, reveal the dynamical trajectories of de-excitation via a conical intersection, yield the key parameters of the conical intersection controlling the de-excitation process, and elucidate the topography of the electronic potential energy surfaces involved.

Molecular hydrogen (H2) is considered an eco-friendly future energy-carrier and an alternative to fossil fuel1 and thus, major efforts are directed towards identifying efficient and economical hydrogen catalysts.2,3 Efficient hydrogen catalysis is used by many microorganisms, some of them producing H2 from organic materials and others consuming it.4-6 To metabolize H2, these microorganisms use enzymes called hydrogenases.7,8 For the future development of efficient catalysts a detailed analysis of the catalytic mechanisms of such hydrogenases is required and existing analytical techniques could not provide a full understanding.9 Consequently, new analytical technologies are of utmost importance to unravel natures blueprints for highly efficient hydrogen catalysts. Here, we introduce signal-enhanced or hyperpolarized, nuclear magnetic resonance (NMR) to study hydrogenases under turnover conditions. So far undiscovered hydrogen species of the catalytic cycle of [Fe]-hydrogenases, are revealed and thus, extend the knowledge regarding this class of enzymes. These findings pave new pathways for the exploration of novel hydrogen metabolisms in vivo. We furthermore envision that the results contribute to the rational design of future catalysts to solve energy challenges of our society.

Sickle cell disease (SCD) is an autosomal recessive genetic disease caused by the Glu6Val mutation in the {beta} chain (Hb) of the oxygen-carrying hemoglobin protein in sicklemia patients. In the molecular pathogenesis of SCD, the sickle hemoglobin (Hb-S) polymerization is a major driver for structural deformation of red blood cells, i.e. red blood cell (RBC) sickling. Biophysically, it still remains elusive how this SCD-linked E6V mutation leads to Hb-S polymerization in RBC sickling. Therefore, with a comprehensive set of analysis of experimental Hb structures, this letter highlights electrostatic repulsion as a key biophysical mechanism of Hb-S polymerization in RBC sickling, which provides atomic-level insights into the functional impact of the SCD-linked E6V substitution from a biophysical point of view.\n\nSIGNIFICANCEDuring the past 25 years, a total of 104 Hb-related structures have been deposited in PDB. For the first time, this article presents a comprehensive set of electrostatic analysis of the 104 experimental structures, highlighting electrostatic repulsion as a fundamental biophysical mechanism for Hb-S polymerization in RBC sickling. The structural and electrostatic analysis here also provides biophysical insights into the functional impact of the SCD-linked E6V substitution.

The thermodynamics of molecular recognition by proteins is a central determinant of complex biochemistry. For over a half-century detailed cryogenic structures have provided deep insight into the energetic contributions to ligand binding by proteins1. More recently, a dynamical proxy based on NMR-relaxation methods has revealed an unexpected richness in the contributions of conformational entropy to the thermodynamics of ligand binding2,3,4,5. There remains, however, a discomforting absence of an understanding of the structural origins of fast internal motion and the conformational entropy that this motion represents. Here we report the pressure-dependence of fast internal motion within the ribonuclease barnase and its complex with the protein barstar. Distinctive clustering of the pressure sensitivity correlates with the presence of small packing defects or voids surrounding affected side chains. Prompted by this observation, we performed an analysis of the voids surrounding over 2,500 methyl-bearing side chains having experimentally determined order parameters. We find that changes in unoccupied volume as small as a single water molecule surrounding buried side chains greatly affects motion on the subnanosecond timescale. The discovered relationship begins to permit construction of a united view of the relationship between changes in the internal energy, as exposed by detailed structural analysis, and the conformational entropy, as represented by fast internal motion, in the thermodynamics of protein function.

The recognition of carbohydrates by lectins play key roles in diverse cellular processes such as cellular adhesion, proliferation and apoptosis which makes it a promising therapeutic target against cancers. One of the most functionally active lectins, galectin-3 is distinctively known for its specific binding affinity towards {beta}-galactoside. Despite the prevalence of high-resolution crystallographic structures, the mechanistic basis and the molecular determinants of the sugar recognition process by galectin-3 are currently elusive. Here we address this question by capturing the complete dynamical binding process of human galectin-3 with its native ligand N-acetyllactosamine (LacNAc) and one of its synthetic derivatives by unbiased Molecular Dynamics simulation. In our simulations, both the natural ligand LacNAc and its synthetic derivative, initially solvated in water, diffuse around the protein and eventually recognise the designated binding site at the S-side of galectin-3, in crystallographic precision and identifies key metastable intermediate ligand-states around the galectin on their course to eventual binding. The simulations highlight that the origin of the experimentally observed multi-fold efficacy of synthetically designed ligand-derivative over its native natural ligand LacNAc lies in the derivatives relatively longer residence time in the bound pocket. A kinetic analysis demonstrates that the LacNAc-derivative would be more resilient compared to the parent ligand against unbinding from the protein binding site. In particular, the analysis identifies that interactions of the binding pocket residues Trp181, Arg144 and Arg162 with the tetrafuorophenyl ring of the derivative as the key determinant for the synthetic ligand to latch into the pocket. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=149 SRC="FIGDIR/small/449218v1_ufig1.gif" ALT="Figure 1"> View larger version (58K): org.highwire.dtl.DTLVardef@1fd3daaorg.highwire.dtl.DTLVardef@db8efaorg.highwire.dtl.DTLVardef@87cacforg.highwire.dtl.DTLVardef@113f0d2_HPS_FORMAT_FIGEXP M_FIG C_FIG

Intermolecular spin relaxation by translational motion of spin pairs have been widely used to study properties of the biomolecules in liquids. Notably, solvent paramagnetic relaxation enhancement (sPRE) arising from paramagnetic cosolutes has gained attentions for various applications, including the structural refinement of intrinsically disordered proteins, cosolute-induced protein denaturation, and the characterization of residue-specific effective near-surface electrostatic potentials (ENS). Among these applications, the transverse sPRE rate known as {Gamma} 2 has been predominantly been interpreted empirically as being proportional to <r-6>norm. In this study, we present a rigorous theoretical interpretation of {Gamma} 2 that it is instead proportional to <r-4>norm and provide explicit formula for calculating <r-4>norm without any adjustable parameters. This interpretation is independent of the type or strength of interactions and can be broadly applied, including to the precise interpretation of ENS.

Understanding weak interactions in protein-ligand complexes is essential for advancing drug design. Here, we combine experimental and quantum mechanical approaches to study the streptavidin-biotin complex, one of the strongest known protein-ligand binders. Using a monomeric streptavidin mutant, we analyze 1H NMR chemical shift perturbations (CSPs) of biotin upon binding, identifying unprecedented upfield shifts of up to -3.2 ppm. Quantum chemical calculations attribute these shifts primarily to aromatic ring currents, with additional contributions from charge transfer effects linked to weak interactions. The agreement between experimental and computed chemical shifts validated the X-ray structure as a reliable basis for detailed computational analyses. Energy decomposition analysis reveals that electrostatics dominate the biotin-streptavidin interaction, complemented by significant orbital and dispersion contributions. Notably, weak non-covalent interactions--such as CH{middle dot} {middle dot} {middle dot} S, CH{middle dot} {middle dot} {middle dot}{pi} , and CH{middle dot} {middle dot} {middle dot} HC contacts--driven by London dispersion forces, contribute [~]44% to the complexs stability.

It has been an established idea in recent years that protein is a physiochemically connected network. Allostery, understood in this new context, is a manifestation of residue communicating between remote sites in this network, and hence a rising interest to identify functionally relevant communication pathways and the frequent communicators within. However, there have been limited computationally trackable general methods to discover proteins allosteric sites in atomistic resolution with good accuracy. In this study, we devised a time-dependent linear response theory (td-LRT) integrating intrinsic protein dynamics and perturbation forces that excite proteins temporary reconfiguration at the non-equilibrium state, to describe atom-specific time responses as the propagating mechanical signals and discover that the most frequent remote communicators can be important allosteric sites, mutation of which could deteriorate the hydride transfer rate in DHFR by 3 orders. The preferred directionality of the signal propagation can be inferred from the asymmetric connection matrix (CM), where the coupling strength between a pair of residues is suggested by their communication score (CS) in the CM, which is found consistent with experimentally characterized nonadditivity of double mutants. Also, the intramolecular communication centers (ICCs), having high CSs, are found evolutionarily conserved, suggesting their biological importance. We also identify spatially clustered top ICCs as the newly found allosteric site in ATG4B. Among 2016 FDA-approved drugs screened to target the site, two interacting with the site most favorably, confirmed by MD simulations, are found to inhibit ATG4B biochemically and be tumor suppressive in colorectal, pancreatic and breast cancer cell lines with an observed additive therapeutic effect when co-used with an active-site inhibitor.

Aberrant misfolding and progressive aggregation of the intrinsically disordered protein (IDP), -synuclein, are associated with the etiology of several neurodegenerative diseases. However, the structurally heterogeneous ensemble of this IDP and lack of a well-defined binding pocket make it difficult to probe the druggability of -synuclein. Here, by building a comprehensive statistical model of the fuzzy ensemble of a millisecond-long atomistic simulation trajectory of monomeric -synuclein interacting with the small-molecule drug fasudil, we identify exhaustive sets of metastable binding-competent states of -synuclein. The model reveals that the interaction with the drug primes this IDP to explore both more compact and more extended conformational sub-ensemble than those in neat water, thereby broadening its structural repertoire in presence of small-molecule via an entropy expansion mechanism. Subsequent simulation of the dimerisation process shows that similar motif of entropic-expansion mechanism helps fasudil to retard the self-aggregation propensity of -synuclein via trapping it into multiple distinct states of diverse compaction featuring aggregation-resistant long-range interactions. Furthermore, small-molecule binding interactions in dimerisation-competent relatively extended states have a screening effect that hinders the formation of stable dimer contacts. Together, the investigation demonstrates the ability of small-molecules to have an ensemble-modulatory effect on IDPs that can be effectively utilised in therapeutic strategies probing aggregation-related diseases.

Cyanobacteria carry out photosynthetic light-energy conversion using phycobiliproteins for light harvesting and the chlorophyll-rich photosystems for photochemistry. While most cyanobacteria only absorb visible photons, some of them can acclimate to harvest far-red light (FRL, 700-800 nm) by integrating chlorophyll f and d in their photosystems and producing red-shifted allophycocyanin. Chlorophyll f insertion enables the photosystems to use FRL but slows down charge separation, reducing photosynthetic efficiency. Here we demonstrate with time-resolved fluorescence spectroscopy that charge separation in chlorophyll-f-containing Photosystem II becomes faster in the presence of red-shifted allophycocyanin antennas. This is different from all known photosynthetic systems, where additional light-harvesting complexes slow down charge separation. Based on the available structural information, we propose a model for the connectivity between the phycobiliproteins and Photosystem II that qualitatively accounts for our spectroscopic data. This unique design is probably important for these cyanobacteria to efficiently switch between visible and far-red light.

Apolipoprotein E (ApoE), a major determinant protein for lipid-metabolism, actively participates in lipid transport in central nervous system via high-affinity interaction with lipoprotein receptor LDLR. Prior evidences indicate that the phospholipids first need to assemble around apoE, before the protein can recognise its receptor. However, despite multiple attempts via spectroscopic and biochemical investigations, it is unclear what are the impact of lipid assembly on globular structure of apoE. Here, using a combination of all-atom and coarse-grained molecular dynamics simulations, we demonstrate that, an otherwise compact tertiary fold of monomeric apoE3 spontaneously unwraps in an aqueous phospholipid solution in two distinct stages. Interestingly, these structural reorganizations are triggered by an initial localised binding of lipid molecules to the C-terminal domain of the protein, which induce a rapid separation of C-terminal domain of apoE3 from the rest of its tertiary fold. This is followed by a slow lipid-induced inter-helix separation event within the N-terminal domain of the protein, as seen in an extensively long coarse-grained simulation. Remarkably, the resultant complex takes the shape of an open conformation of lipid-stabilised unwrapped protein, which intriguingly coincides with an earlier proposal by small-angle X-ray scattering (SAXS) experiment. The lipid-binding activity and the lipid-induced protein conformation are found to be robust across a monomeric mutant and wild-type sequence of apoE3. The open complex derived in coarse-grained simulation retains its structural morphology after reverse-mapping to all-atom representation. Collectively, the investigation puts forward a plausible structure of currently elusive conformationally activated state of apoE3, which is primed for recognition by lipoprotein receptor and can be exploited for eventual lipid transport. Table of Content O_FIG O_LINKSMALLFIG WIDTH=199 HEIGHT=200 SRC="FIGDIR/small/255299v2_ufig1.gif" ALT="Figure 1"> View larger version (67K): org.highwire.dtl.DTLVardef@1551278org.highwire.dtl.DTLVardef@2f85fdorg.highwire.dtl.DTLVardef@5c90aorg.highwire.dtl.DTLVardef@1d1ddb2_HPS_FORMAT_FIGEXP M_FIG C_FIG

An increasing number of ligand-bound membrane protein structures reveal ligand-binding sites on the lipid-exposed surface of the protein within the membrane bilayer. Binding events to such sites have previously been studied using molecular dynamics (MD) simulations and experiments in cases such as calcium-gated potassium channels1 and sodium channels2. The proposed binding mechanism is that these ligands partition into the membrane to gain access to their binding site. What is currently unavailable is what the thermodynamic and kinetic contributions of the ligand-membrane and ligand-protein interactions are to the overall binding event. Here, we used MD simulations and enhanced sampling methods to study the membrane partitioning of a DHP calcium channel antagonist, nifedipine, from the voltage-gated calcium channel Cav1.1. We present that drug-membrane interactions occur on a much faster timescale than the overall binding of nifedipine to Cav1.1.